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Abstract 

Background: There is strong but mostly circumstantial evidence that genetic factors modulate the severity of 
influenza infection in humans. Using genetically diverse but fully inbred strains of mice it has been shown that 
host sequence variants have a strong influence on the severity of influenza A disease progression. In particular, 
C57BL/6J, the most widely used mouse strain in biomedical research, is comparatively resistant. In contrast, 
DBA/2J is highly susceptible. 

Results: To map regions of the genome responsible for differences in influenza susceptibility, we infected a 
family of 53 BXD-type lines derived from a cross between C57BL/6J and DBA/2J strains with influenza A virus 
(PR8, HlNl). We monitored body weight, survival, and mean time to death for 13 days after infection. Qivr5 
(quantitative trait for influenza virus resistance on chromosome 5) was the largest and most significant QTL for 
weight loss. The effect of Qivr5 was detectable on day 2 post infection, but was most pronounced on days 5 and 6. 
Survival rate mapped to Qivr5, but additionally revealed a second significant locus on chromosome 19 {Qivrl9). 
Analysis of mean time to death affirmed both Qivr5 and O/Vr/9. In addition, we observed several regions of the 
genome with suggestive linkage. There are potentially complex combinatorial interactions of the parental alleles 
among loci. Analysis of multiple gene expression data sets and sequence variants in these strains highlights about 
30 strong candidate genes across all loci that may control influenza A susceptibility and resistance. 

Conclusions: We have mapped influenza susceptibility loci to chromosomes 2, 5, 16, 17, and 19. Body weight 
and survival loci have a time-dependent profile that presumably reflects the temporal dynamic of the response to 
infection. We highlight candidate genes in the respective intervals and review their possible biological function 
during infection. 



Background 

Influenza A virus represents a major health threat to 
humans. The 1918 HlNl pandemic caused at about 30 
to 50 million deaths [1]. Seasonal influenza epidemics 
cause high economic loss, morbidity and deaths every 
year [2]. The course and outcome of an influenza A 
virus infection is influenced by viral and host factors. 
Host risk factors, like obesity or pregnancy, became evi- 
dent during the recent swine flu pandemics [3,4]. 
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Furthermore, genetic factors in humans associated with 
a higher susceptibility to influenza infections and severe 
disease outcome have been suspected for the 1918 
pandemics, as well as the H5N1 human infections [5-7]. 
Recently, the importance of IFITM3 as a crucial factor 
for host susceptibility has been demonstrated in mice 
and humans [8]. 

The importance of host factors to host susceptibility 
and resistance has been demonstrated clearly in animal 
models. We and others have shown in mouse infection 
models that the susceptibility of the host to influenza A 
infection strongly depends on the genetic background 
[9-17]. In particular, DBA/2J mice are highly susceptible 
to many influenza A virus subtypes, including those that 
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were directly derived from human isolates without prior 
adaptation to the mouse [9,13,16-18]. In contrast, 
C57BL/6J mice are more resistant. After infection with 
mouse-adapted HlNl (PR8M virus), DBA/2J mice loss 
weight very rapidly and die within 5-7 days post infec- 
tion (p.i.)> whereas C57BL/6J mice loss weight until days 
6-8 after infection and regain their initial weight by 
14 days p.i. [9,18]. Viral load in the lungs of DBA/2J 
infected mice is much higher and lung pathology is very 
severe compared to infected C57BL/6J mice. Also, the 
production of chemokines and cytokines is much higher 
in DBA/2J mice [9,18]. 

However, the genomic regions that are responsible for 
the differential response after infection with HlNl have 
not been determined. Therefore, we used a large family 
of BXD type recombinant inbred strains generated by 
crossing C57BL/6J (resistant) to DBA/2J (susceptible) to 
map genetic loci that modulate disease severity. The 
BXD genetic reference population (GRP) is made up of a 
set of progeny strains, each with a defined and fixed gen- 
etic architecture. It is one of the largest families of 
strains, consisting of about 80 fully inbred strains [19,20] 
available from the Jackson Laboratory and a new set of 
80 additional lines that are still in production at the Uni- 
versity of Tennessee. Individuals within each single 
strain are essentially isogenic (except for the sex chro- 
mosomes) and genotypes for the entire family, including 
most of the new strains, are known [21]. Genetic vari- 
ation among this family has be exploited extensively in 
the past to systematically study the genetics of many 
traits (for examples of phenotypes see the GeneNetwork 
database [22]). 

Here, we infected over 50 of the BXD strains with 
influenza A HlNl virus and monitored body weight, 
survival, and mean time to death for the following 
13 days post infection. We identified two significant and 
several suggestive loci peaks for all three traits. All 
showed a time-dependent appearance. Data mining of 
the intervals revealed several candidate genes, several 
of which may be important for the host response to 
influenza A virus infection. 

Results 

Susceptibility to influenza A virus after experimental 
infection of BXD mouse strains is highly variable 

We infected 53 recombinant inbred strains of the BXD 
population plus the parental strains C57BL/6J, DBA/2J 
and B6D2F1 mice with mouse-adapted HlNl virus 
(PR8M, HlNl [18]) and followed body weight and 
survival over the next 13 days. Mice that lost more than 
25% of their starting weight were sacrificed and also 
recorded as dead. After infection, body weight changes 
were highly variable between the BXD strains over 
the period of 13 days post infection (Figure 1, Additional 



file 1: Table SI). Similarly, survival rates (Additional 
file 2: Figure SIA, B) and mean time to death (MTTD; 
Additional file 2: Figure SIC) were highly variable 
between the different BXD strains. Furthermore, the 
trait percent survival showed a strain-dependent progres- 
sion over time. 

For body weight loss and survival, three different 
phenotypic response groups can be defined (Figure 1). In 
the first (Figure lA), all infected mice within a strain sur- 
vived, in the second (Figure IB) a majority but not all 
individuals within a strain survived, and in the third 
group (Figure IC), a majority died. Most remarkably, in 
the first group four strains— BXD9, BXD13, BXD43, 
BXD97 — were highly resistant indicating that the infec- 
tion may not cause any major pathology (Figure lA). In 
contrast, BXD28 belonging to the third group, lost body 
weight much more rapidly than even the highly suscep- 
tible DBA/2 J parent (Figure IC). These results illustrate a 
large variation of responses within the BXD family. 

By day 7 p.i., all infected mice had succumbed to 
infection in 11 strains, whereas 14 others exhibited 
more limited mortality (Additional file 2: Figure SIA). The 
incidence of mortality increased in some strains from day 
8 p.i. through day 11 but not thereafter (Additional file 2: 
Figure SIB). For the MTTD phenotype, 17 strains showed 
no mortality after infection, similar to the resistant 
C57BL/6J parent. Three strains— BXD28, BXD18, 
BXD103— exhibited a MTTD that was even shorter 
than for the susceptible DBA/2J parent (Additional file 2: 
Figure SIC). Although this study was conducted over a 
period of approximately three years, and although mice 
were received from different sources and different experi- 
menters performed the infection experiments, we did not 
note any significant influence of these potential confounds 
and cofactors. 

Principal Component Analysis (PCA) of strain mean 
body weight loss from day 1 until day 7 was carried out 
to reduce the number and redundancy of related mea- 
surements, as well as to evaluate whether this set of mea- 
surements can be broken down into statistically and 
genetically independent processes. This analysis revealed 
one major component, PCI that explains most of the weight 
loss among strains (more than 80%, see Additional file 3: 
Figure S2). This component corresponds to differences 
in loss from day 3 to 7 p.i. (Figure 2). In contrast, PC2 
accounts for only about 10% of the variance in weight loss 
(Additional file 3: Figure S2). This component corresponds 
only to very early weight loss after virus inoculation, on 
days 1 and 2 p.i. (Figure 2), well before any appreciable 
virus replication has occurred. Thus, PC2 is highly likely to 
be a technical, but still interesting effect associated with 
the stressful procedure of intranasal inoculation and full 
anesthesia. PCI most likely reflects the biologically relevant 
variation among strains in response to viral replication. 
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Figure 1 BXD strains exiiibit variable Icinetics of weight loss and survival after infection with Influenza A virus. Mice from 53 BXD and 
parental strains were infected intra-nasally witli 2x10^ FFU of PR8 virus. Weiglit loss and survival of infected mice was followed over a period of 
13 days. Mortality includes mice that were sacrificed because they had lost more than 25% of body weight. Three phenotypic groups can be 
distinguished: in the first group (A), all infected mice within a given BXD strain survived, in the second group (B) less than 50% of infected mice 
within a given BXD strain died, in the third group (C), more than 50% of infected mice within a given BXD strain died. From the weight loss 
curve of the second and third group it also becomes obvious that non-surviving mice were all approaching the 75% body weight loss endpoint 
before dying. 



Analysis of body weight loss revealed a significant QTL 
on chromosome 5 and several suggestive QTLs with 
time-dependent effects 

We mapped body weight loss following infection day by 
day. For purposes of analysis, body weights of mice that 
had died or that were euthanized were assigned a weight 
equal to 75% of their initial weight. It should be noted 



that mice which died continuously lost weight and were 
close to 75% body weight loss before they were found 
dead (Figure IB, C). After day 7, surviving mice started 
to gain weight (Figure lA, B). For this reason, we limited 
our analyses of body weight traits to the period of day 1 
to day 7 p.i. in order to avoid mixing data from dead 
with recovering mice. Interval mapping of body weight 
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Figure 2 Principal component analysis of all body weight traits 
reveals two distinct groups. PCA loading plot of principal 
components PCI and PC2 (Supplement Figure 2) for body weight 
traits form day 1 to day 7 p.i. shows two distinct groups. The 
variation of body weight loss at day 3 to day 7 (ID 13007 to ID 
1301 1) are mainly represented by factor 1. On the other hand, 
variation of body weight loss at day 1 (ID 13005) and 2 (ID 13006) 
are most likely explained by factor 2. 

V J 

loss detected a significant time-dependent locus on 
chromosome 5 that we named 'QTL for influenza virus 
resistance on chromosome 5' {QivrS) adopting the no- 
menclature proposed by [12]. This QTL exhibited 
genome-wide significance on days 5 (LRS: 19.0, effect 
size: 30%) and 6 (LRS: 19.4, effect size: 31%, Figure 3). 
However, the effect of QivrS, although weaker, can also 
be detected on days 2, 3, and 4 and 7 (Figure 3). Most 
interestingly; the resistance allele at QivrS is inherited 
from the nominally sensitive DBA/2J parental strain, 
illustrating the genetic complexity of the influenza 
response. Suggestive loci (LRS values between 10 and 
15) map to chromosomes 2, 6, 9, proximal and distal 16, 
and chromosome 17. The effect of these loci was also 
time-dependent (Figure 3A-G). 

Another suggestive QTL peak was found on chromo- 
some 10 on day 1 (LRS: 14.1, effect size: 23%, 
Figure 3A). Its effect is lower at day 2 and not apparent 
at later days. These observations are in accordance with 
the PCA (Figure 2) that reveals two separate time- 
dependent influences on body weight variance among 
strains. The QTL appeared at an early time point 
after infection, when virus replication has just begun 
and strong inflammatory host responses are not yet 
evident [9,18,23]. These observations indicate that the 
effect is most likely related to the experimental proto- 
col, namely the stress to anesthesia and intra-nasal 
application as well as treatment recovery. Treatment- 
dependent QTLs were described previously [24]. It is 
worth to note that mock-infection of the parental 
DBA/2J and C57BL/6J mice did not lead to a lasting 
body weight loss over a longer time interval except for 



a slight drop in body weight on day 1 p.i. (Additional 
file 4: Figure S3). 

We also performed a QTL analysis for the PCI and 
PC2 described above. PCI detects a significant QTL on 
chromosome 5 (LRS: 19, effect size: 30%, Figure 4A) as 
well as suggestive QTL peaks on proximal chromosome 
16 and on chromosome 17 (Figure 4A). Thus, the PCI 
confirmed the significant QTL on chromosome 5 found 
with the daily body weight loss traits. Interval mapping 
of PC2 detected no significant QTL and one suggestive 
QTL on chromosome 13 (Figure 4B) which was not seen 
in any other trait. 

Survival rate and mean time to death traits confirmed 
the QTL on chromosome 5 and detected another 
significant QTL on chromosome 19 

The analysis of survival rate traits on days 7, 8 and 11 p.i. 
revealed significant peaks on chromosomes 5 and 19 and 
suggestive peaks on chromosomes 1, 2, 3, 10, 16 and 17 
(Figure 5A-C). The effects of these QTLs were again 
time-dependent. QivrS was significant at day 7 (LRS: 
24.8, effect size: 37%) and day 8 p.i. (LRS: 19.7, 
effect size: 31% Figure 5 A, B) and its effect was still evi- 
dent from day 9 until 11 p.i. (Additional file 5: Figure S4 
and Figure 5C). Similarly, another significant QTL with 
time-dependent effects was observed on chromosome 19 
{Qivrl9) at day 8 (LRS: 19.4; effect size: 29%) that could 
also be detected at day 7 and days 9-11 (Figure 5A-C 
and Additional file 5: Figure S4). QivrS and Qivrl9 
resulted from a positive influence of DBA/2J on survival, 
whereas Qivrl6 (distal locus) and Qivrl7-2 resulted from 
a positive influence of C57BL/6J alleles to increase sur- 
vival. Furthermore, MTTD analysis confirmed QivrS as a 
significant QTL (LRS: 20.5, effect size: 32%), and Qivrl9 
was almost significant for this trait (Figure 5D). In 
addition, suggestive QTLs were found on chromosomes 
2 and 17 for the MTTD trait. 

In conclusion, survival and MTTD traits confirmed the 
significant QTL on chromosome 5, revealed an additional 
significant QTL on chromosome 19 and several suggest- 
ive QTLs. All QTLs showed a time-dependent effect. 

Composite interval and pair-scan mapping indicates 
various interactions of QTLs 

The strongest QTLs map to chromosomes 5 and 19. We 
therefore performed composite interval mapping in 
which the contributions of these strong QTLs were fac- 
tored out to reveal possible secondary QTLs. When we 
controlled for Qivrl9, the linkage to distal chromosome 
16 increased and became nearly significant for both sur- 
vival at day 8 and MTTD (Figure 6A, B). The linkage to 
chromosome 2 increased slightly. However, control for 
the contribution of QivrS did not reveal any additional 
loci linked to survival or MTTD traits (Figure 6C, D). 
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Figure 3 Genome-wide linlcage analysis for body weight loss after H1N1 infection identifies time-dependent QTLs. Analysis of body 
weight from day 1 (A), day 2 (B), day 3 (C), day 4 (D), day 5 (E), day 6 (F) and day 7 (G) after PR8 infection revealed time-dependent QTLs. 
A significant QTL is located on chromosome 5 and suggestive QTLs are observed on chromosomes 2, 6, 9, 10, 16 and 17. Dead or euthanized 
mice were included for the analysis by manually setting their body weight to 75%. Interval mapping is shown across the entire genome as blue 
line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows 
mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines 
mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line). 



Finally, we analyzed the joint additive effects of QTLs 
described above as well as possible epistatic interactions 
among these QTLs. When using a full two-locus model 
(Trait Variance = Ql + Q2 + Ql x Q2), we found two 
potential interactions between QivrS and a new locus 
on chromosome 9 as well as between QivrS and Qivrl9, 
respectively, for the survival trait at day 8 (Figure 7A). 
The inclusion of an interaction term increases the LRS 
by 23 (Additional file 6: Figure S5). For both QivrS - 
chromosome 9 interaction (Figure 7B) and QivrS - 
Qivrl9 interaction (Figure 7C), the allele combination 
DBA/2J / DBA/2J showed highest and the combination 
C57BL/6J / C57BL/6J lowest survival scores. All differ- 
ences were highly significant. 



Analysis of QTL regions identified several candidate 
genes that may contribute to host susceptibility or 
resistance 

In total, the mapping studies revealed five QTLs on 
chromosomes 2, 5, distal 16, 17 and 19 {Qivr2'2, QivrS, 
Qivrl6, Qivrl7'2 and Qivrl9) that did merit further 
analysis because they were consistently observed in at 
least two traits and exerted an effect on at least two dif- 
ferent days p.i. 

For the candidate gene searches, we defined the critical 
regions of the QTLs manually by considering peak 
height, its shoulder and bootstraps as guidance (Figure 8). 
We then investigated these Qivr intervals for potential 
candidate genes that may be causal for the studied traits 
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Figure 4 Genome-wide linlcage scan for PCI and PC2 body weigiit loss. Interval mapping of principle component PCI (A) for body weight 
loss traits at day 1 to 7 after PR8 infection shows significant QTLs on chromosome 5 and suggestive QTLs on chromosomes 16 and 17. Interval 
mapping of PC2 (B) for body weight loss traits reveals no significant QTLs and a suggestive QTL on chromosome 13. Interval mapping is shown 
across the entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BI76J alleles increase 
trait values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS 
score. The horizontal lines mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line). 



(Additional file 7: Figure S6). We first used the QTLminer 
tool in GeneNetwork [25] to identify all genes in the QTL 
intervals and to obtain associated GO terms. Next, we 
select all annotated genes within the QTL intervals that 
were expressed between day 1 and day 60 after infection 
of C57BL/6J mice with PR8M virus (the latter data set is 
derived from a separate study (Pommerenke et al, PLoS 
ONE, in press). Subsequently, these genes were further 
characterized for the following attributes: genes that were 
up- or up-regulated in infected lungs by at least 1.5-fold 
in C57BL/6J after PR8M infection, genes carrying an in- 
sertion or deletion or a non-synonymous nucleotide 
change in the open reading frame, genes exhibiting a cis- 
expression QTL in the lung [26], and genes that were dif- 
ferentially expressed by at least 1.5-fold in C57BL/6J and 
DBA/2J mice in infected lungs between day 1 and day 8 



p.i. [23] and Pommerenke et al, PLoS ONE, in press). 
The strategy of the QTL mining is shown as flow chart in 
the Additional file 7: Figure S6, and the final results are 
presented in Table 1. Furthermore, the attributes of all 
genes located in the QTL regions are listed in detail in 
the Additional file 8: Table S2). Using these combined 
attributes, we identified 31 genes as the most likely candi- 
dates to regulate the traits controlled by Qivr2-2, QivrS, 
Qivrl6, Qivrl7-2 and Qivrl9 (listed in Table 2). These 
genes were fiirther evaluated based on their known func- 
tion from the literature and phenotypes in knock-out 
mouse mutants (see discussion). 

Discussion 

DBA/2J and C57BL/6J mice have been shown pre- 
viously to differ largely in their susceptibility to HlNl 
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Figure 5 Genome-wide linlcage scan for survival rate and mean time to death confirms QTLs on chromosomes 5 and 19. Interval 
mapping of survival rates (in %) after PR8 infection at day 7 (A), day 8 (B) and day 1 1 (C) show significant QTLs on chromosomes 5 and 19. 
Suggestive QTLs are observed on chromosomes 2, 3, 10, 16 and 17. Interval mapping using mean time to death confirmed the significant QTL on 
chromosome 5 and revealed suggestive QTLs on chromosomes 2, 17 and 19 (D). Interval mapping is shown across the entire genome as blue 
line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows 
mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines 
mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line). 
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Figure 6 Composite interval mapping indicates interactions of QTLs on chromosomes 5 and 19. The influence of the markers rsl 3483633 
on chromosome 19 (A, B) and rsl 3478587 on chromosome 5 (C, D) was factored out and the residual LRS were calculated. Two different traits, 
survival day 8 (A, C) and mean time of death (B, D) were tested by using the web tool of GeneNetwork Interval mapping is shown across the 
entire genome as blue line representing the LRS scores; green line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait 
values. Upper x-axis shows mouse chromosomes, lower x-axis shows physical map in mega bases for each chromosome, y-axis represents LRS 
score. The horizontal lines mark the genome-wide significant thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line). 
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Figure 7 Interactions are observed between several QTLs. The graph displays pair-scan results for the trait percent survival at day 8 p.i. 
(IDl 3002) (A). The upper left half of the plot highlights any epistatic interactions, the lower right half provides a summary of LRS of the full 
model, representing cumulative effects of linear and non-linear terms based on the model Trait Variance = Q1 -fQ2 + Q1 x Q2 '. Subsequently, 
BXD lines were grouped according to their allele combinations and trait values were compared between groups (B, C). The largest differences are 
indicated by asterisks. They are highly significant between the respective groups (p< 0.0001). For the various loci, trait values for the following 
markers were used: Q/\/r5: rsl 3478587, chromosome 9: rs6191976, Q/Vr/9: rsl 3483633. Number of lines and p-values for Q/Vr5 - chromosome 9: 
D2/D2: n = 1 7, D2/B6: n = 1 7, B6/D2: n = 6, B6/B6 n = 1 2; D2/D2 vs B6/B6: p < 0.0001 , D2/B6 vs B6/B6: p = 0.003. Number of lines and p-values for 
QmS - Oivrl9: D2/D2: n = 18, D2/B6: n = 16, B6/D2: n=4, B6/B6: n = 14; p-values: D2/D2 vs B6/B6: p< 0.0001, D2/B6 vs B6/B6: p = 0.0008. 
Abbreviations: D2: DBA/2J, B6: C57BL/6J. 



(PR8M) influenza A virus [9,18]. Here, we expanded 
these studies and utilized the BXD recombinant inbred 
set of mouse strains to map the genomic regions that 
are responsible for differences in these two mouse 
strains. We monitored three phenotypic traits, body 
weight over time, survival over time and mean time to 
death to identify quantitative trait for influenza resist- 
ance. Two significant QTLs, QivrS and Qivrl9, were 
found on chromosomes 5 and 19, respectively. Fur- 
thermore several suggestive QTLs, Qivr2'2, Qivrl6 and 
Qivrl7-2 were observed in at least two traits and at 
two days on chromosomes 2, 16 and 17, respectively. 
Composite mapping revealed an additional almost sig- 
nificant QTL at distal chromosome 16, Qivrl6, 

A similar analysis for host resistance to influenza has 
been performed previously after infecting 66 BXD 
strains with H5N1 influenza virus. This study reported 
three significant QTLs on chromosomes 2, 7, and 17 
[12]. Thus, none of these significant QTLs overlaps with 
the QTLs identified in our analysis. Five of the strains 
that were resistant (all infected mice survived) in our 



study were also resistant in the study of [12] where a 
total of 14 strains were found to be resistant. Five strains 
that were highly susceptible in our study (100% of 
infected mice died) were also highly susceptible in the 
study by [12] of a total of 26 susceptible strains. Further- 
more, five strains that were resistant in our study were 
susceptible in the study by [12]. Thus, there is also 
not much overlap between the two studies with respect 
of susceptible and resistant strains. The differences 
between the two studies are most likely explained by the 
use of two different influenza virus subtypes. The HlNl 
virus from our study represents a subtype with a mono- 
basic hemagglutinin (HA) cleavage site, whereas the 
H5N1 which was used in the study by Boon et al. is a 
subtype with a polybasic HA cleavage site. The cellular 
tropism of these two subtypes for virus replication and 
processing is quite different, because monobasic viruses 
are dependent on cell-specific proteases for the pro- 
cessing of the HA whereas polybasic subtypes can be 
processed by more ubiquitously expressed host pro- 
teases, e.g. [27-32]. Therefore, the contribution of host 
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(See figure on previous page.) 

Figure 8 Detailed maps of QTLs contributing to host susceptibility and resistance to PR8 influenza A infection. Interval mapping was 
performed for chromosomes 2 (A), 5 (B), 16 (C), 17 (D) and 19 (E). The critical intervals were selected based on the peak shapes and the 
bootstrap signals: Oivr2-l (56-68 Mb), QivrS (140-153 Mb), Q/Vr/6 (64-78 Mb), Oivrl7-2 (30-44 Mb) and Q/Vr/9 (37-45 Mb). The Qivr intervals are 
delineated by red arrows. The yellow bars represent the frequencies of peak LRS using bootstrap analysis. The multicolored chatters along the 
top of the graph are hyperlinks to sites with additional genetic and sequence information. The orange chatter along the x-axis indicates the 
density of SNPs present in the BXD strains. Interval mapping is shown across the entire genome as blue line representing the LRS scores; green 
line: DBA/2J alleles increase trait values; red line: C57BL/6J alleles increase trait values. Upper x-axis shows mouse chromosomes, lower x-axis 
shows physical map in mega bases for each chromosome, y-axis represents LRS score. The horizontal lines mark the genome-wide significant 
thresholds at p<0.05 (red line) and suggestive thresholds at p<0.37 (gray line). 



factors to susceptibility may be different between HI 
and H5 containing virus subtypes. 

Another study described the genetic mapping of sus- 
ceptibiUty and resistance factors after infecting a panel 
of 29 AxB / BxA congenic strains with a mouse-adapted 
H3N2 influenza virus [33]. The AxB / BxA congenic 
strains were generated from a cross of susceptible A/J 
and resistant C57BL/6J parental mouse strains. The 
authors found three major QTLs on chromosomes 2, 6 
and 17. The QTL on chromosome 17 overlaps with the 
Qivr 17-2 locus which we found in our study. Further- 
more, the candidate gene Pla2g7 that was identified in 
their study was also detected as candidate gene in our 
analysis (see below). 

The influence of genetic factors determining the host 
response to HlNl influenza virus infections was also 
examined in mice of the pre-Collaborative Cross collec- 
tion [34]. In this study, gene expression levels in extreme 
responders were used to identify expression QTLs 
(eQTL). One gene that exhibited a cis-eQTL, Sikl (salt 
inducible kinase 1), was located in the Qivr 17-2 interval 
from our study, and we also identified it as potential 
quantitative trait gene (Table 2). This gene is associated 
with the GO terms negative regulation of transcription 
from RNA polymerase II promoter, regulation of cell 
differentiation, and protein kinase cascade'. However, no 
specific infection-related fiinctions have been yet described 
for this gene. 

One of the most interesting findings in our study was 
the time-dependent effect of QTLs which we observed 
in the body weight and survival traits. The peak QTLs 
for the two significant QTLs, QivrS and Qivr 19, were 
found at different times p.i., day 6 and day 8, respect- 
ively. In addition, the effects of both QTLs were not only 
evident at the times p.i. where they exerted the signifi- 
cant peak QTL signals but also several days before and 
after the peak. Furthermore, for the suggestive QTLs, 
also time-dependent effects were observed. These results 
suggest that the causal genes underlying different QTLs 
act at different time points of the host defense. 

Most interestingly, QivrS as well as Qivr 19 represent a 
positive influence on body weight, survival and MTTD 
from the DBA/2J haplotype, the susceptible strain. 



These findings indicate that genomic regions from the 
susceptible parent are able to increase resistance when 
combined with the resistant parental genome. We are 
now analyzing several BXD strains that were more 
resistant than the parental C57BL/6J mice in more 
detail. One possible mechanism to explain such an effect 
may be that an activator (secreted ligand or transcription 
factor) is expressed in susceptible DBA/2J mice but the 
corresponding target (receptor or regulated gene) is 
mutated. On the other hand, in C57BL/6J mice, the tar- 
get but not the activator may be mutated. If the wild 
type alleles are now coming together in a BXD strain, 
the functional activator finds its functional target and 
thereby an increased resistance state is achieved. 

Both composite and interaction mapping revealed 
many genetic interactions between C57BL/6J and 
DBA/2 J alleles. Thus, many genomic regions from the 
parental strains are able to contribute to the host re- 
sponse and this effect depends strongly on the allele 
combinations in the respective QTLs. These observa- 
tions may be studied further in double congenic 
mouse lines. 

We subsequently analyzed the five QTL intervals, 
Qivr2-2, QivrS, Qivr 16, Qivr 17-2 and Qivr 19 in more de- 
tail to identify genes that may be causal for resistance or sus- 
ceptibility. In total, 830 genes are located in these intervals. We 
narrowed down the total number of genes to 31 candidates 
(Table 2) by using additional information, such as temporal ex- 
pression after PRM8 infection (Pommerenke et al, PLoS 
ONE, in press), cis-eQTLs in non-infected lungs [26], dif- 
ferences in expression between DBA/2J and C57BL/6J 
[23], and sequence variants in the coding regions. 

QivrS contains the candidate gene Eif2akl (eukaryotic 
translation initiation factor 2 alpha kinase 1) that is a 
member of eIF2alpha kinases which have been associated 
with anti- viral host responses [35]. Boon et al. described 
another eIF2alpha kinase, Eip,ak2 I Pkr (eukaryotic trans- 
lation initiation factor 2-alpha kinase 2), in the Qirvl7 
locus after infection with influenza H5N1 [12]. Eif2ak2 
plays a critical role in modulating immunoglobulin ex- 
pression during RSV infection. In addition Eif2ak2 knock- 
out mouse mutants are more susceptible to influenza 
infections [36,37]. We have initiated the generation of a 
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Table 1 Candidate genes in mapped QTL intervals 



Qivr 


Interval 


Allele No of 
increasing genes in 
survival interval 


No of genes 
expressed 
during infection 


No of genes 
with Indels 
(FS in coding region) 


No of genes with Cis-eQTLs in 
SNPs (non-synonymous non-infected lung 
codons / stop codons) 


Differentially 
expressed 
btw B6/D2 










(up-/down- regulated) 










2-2 


56-68 Mb 


D2 


72 


24 (19) 




4/0 


6 


3 


5 


140-153 Mb 


D2 


179 


83 (57) 




9/0 


13 


10 


16 


64-78 Mb 


B6 


52 


12(7) 


1^ 


5/0 


4 


4 


17-2 


30-44 Mb 


B6 


370 


158 (104) 


1^ 


58 / 4'"^ 


32 


40 


19 


37- 45 Mb 


D2 


161 


51 (30) 


19 


14/ 1^ 


11 


4 



Up- or down-regulated genes were defined as genes tliat exiiibited at least a 1.5 difference change in expression levels in lungs of infected mice compared to the 
non-infected controls. Genes were defined as differentially expressed exhibiting at least a 1.5 difference in expression changes in infected lungs of C57BL/6J and 
DBA/2J. Genes with sequence variations: ^Robol, ^Lstl, ^Brd2, '^H2-Ab1, ^H2-BI, ^Rpp21, ^TctnS, ^Dhdpsl. Abbreviations: FS: frame shift in coding region, SNP: single 
nucleotide polymorphism, cis-eQTL: cis regulated expression QTL, B6: C57BL/6J, D2: DBA/2J. 



congenic mouse lines for the chromosome 5 interval to 
verify and further characterize the effect of this region for 
resistance to influenza infection. 

Qivr2-2 contains two candidate genes, Itgb6 (integrin 
beta 6) and Ifihl (interferon induced with helicase C do- 
main 1), with known functions in the host defense 
to viral infections. Itgb6 mouse knock-out mutants 
exhibit severe pneumonia and an increase in granulocyte 
recruitment to the lung [38]. The protease-activated 
receptor 1 -mediated enhancement of /^^Z?6-dependent 
TGF-beta activation has been proposed to represent one 
mechanism by which activation of the coagulation cas- 
cade contributes to the development of acute lung injury 
[39]. The Ifihl gene is also known as MDA5 (Melanoma 
Differentiation- Associated protein 5). IFIHl is part of 
the RIG-I-like receptor (RLR) family, which function as 
pattern recognition receptors and are activated upon 
binding of virus dsRNA [40]. IFIHl functions as cyto- 
solic receptor that leads to the selective activation of type 
I IFN genes and is indispensable for sustained expression 
of IFN in response to paramyxovirus infection [41,42]. 
Ifihl mutant knock-out mice exhibit an impaired response 
to different viral pathogens [43,44]. 

Qivr 16 contains two potential genes with known 
functions in the host defense and lung function, Robol 
(roundabout homolog 1 (Drosophila)) and Nripl (nuclear 
receptor interacting protein 1). DBA/2 J mice carry a frame 
shift mutation in the Robol gene which might lead to an 
impaired function of the encoded protein. Robol has been 
described to be involved in guidance and migration of axons, 
myoblasts, and leukocytes in vertebrates (e.g. [45-47]J but is 
also expressed in the developing lung [48]. Robol knock- 
out mutants exhibit a delayed lung maturation and bron- 
chial hyperplasia. The latter results suggest that Robol 
may be involved in maintaining proper lung function and 
it may become essential when lung epithelium is 
destroyed during an influenza infection. Nripl/RipMO 
functions as a co-activator for cytokine gene promoter ac- 
tivity via direct protein-protein interactions with the 



NFkappaB subunit ReLA and histone acetylase cAMP- 
responsive element binding protein (CREB) -binding pro- 
tein (CBP) [49]. It is involved in modulating pro- 
inflammatory responses in macrophages [50]. 

Qivr 17-2 represents a positive influence of the C57BL/6 J 
genotype on body weight, survival and MTTD. This 
QTL is located in a gene-rich region which carries many 
genes that are involved in the host immune response, 
in particular the H2 histocompatibility genes which are 
involved in antigen presentation [51]. Therefore, many 
candidate genes are found in the Qivr 17-2 region. The 
Lstl (leukocyte specific transcript 1) gene is of special 
interest because the DBA/2J allele mice carries a single 
nucleotide deletion in the first exon resulting in a frame 
shift of the open reading frame. This mutation most 
likely results in a non-functional Lstl protein in DBA/2 J 
mice. We confirmed the presence of the deletion by se- 
quencing the parental DBA/2J and some BXD strains car- 
rying the DBA/2J allele. The wild type allele was 
confirmed in C57BL/6J mice and in some BXD strains 
carrying the C57BL/6J allele. In humans, LSTl plays a 
role in the regulation of the immune response to in- 
flammatory diseases such as rheumatoid arthritis, mi- 
crobial infection or Rubella vaccine-induced immunity 
[52-55]. Also, Lstl is up-regulated after influenza A in- 
fection in C57BL/6J mice starting at day 2 and exhibits 
a strong peak of expression at day 8 p.i. Pommerenke 
et al, 2012 (Pommerenke, C, E. Wilk, B. Srivastava, A. 
Schulze, N. Novoselova, R. Geffers, and K. Schughart. 
2012. Global transcriptome analysis in influenza-infected 
mouse lungs reveals the kinetics of innate and adaptive 
host immune responses. PLoS ONE. 7:e41169.). Thus, 
the expression profile and known functions of Lstl fit 
well with a possible critical role for the host defense to 
influenza A virus. We initiated the generation of knock- 
out mice to evaluate the role of Lstl in more detail. In 
addition, a second, most interesting candidate, Pla2g7 
(phospholipase A2, group VII (platelet-activating factor acetyl- 
hydrolase, plasma)) was identified in the Qivr 17-2 interval. In 



Nedeiko et al. BMC Genomics 2012, 13:41 1 
http://www.bionnedcentral.conn/1471 -21 64/1 3/41 1 



Page 1 3 of 1 9 



Table 2 GO-terms and functions observed in Icnoclc-out mice of 31 potential candidates from QTL intervals 



Qivr 


Gene 
symbol 


Gene description 


Function 


KO phenotype 


Type 

polymorphism 
in ORF 


Oivr2-2 


Itgb6 


integrin beta 6 


Integrin-mediated signaling 
pathway, inflammatory response, 
cell-matrix adhesion 


Baldness associated with 
macrophage infiltration of skin, 
exaggerated pulmonary 
inflammation, impaired mucosal 
mast cell response to nematode 
infection. 


ns (1) 


Oivr2-2 


Ifitil 


interferon induced witli 
lielicase C domain 1 


Response to virus, innate 
immune response, regulation 
of apoptosis, RIG-l-like receptor 
signaling pathway 


Increased virus-associated morbidity 
and mortality, decreased cytokine 
response to several viral infection. 


ns (5) 


QivrS 


Eif3b 


eul<aryotic translation 
initiation factor 3, 
subunit B 


Translation, translation initiation 


NA 


ns (1) 


OivrS 


Sdkl 


sidel<icl< liomolog 1 
(cliicl<en) 


Cell adhesion 


NA 


ns (1) 


QivrS 


Eif2Qk 1 


eul<aryotic translation 
initiation factor 2 alpha 
kinase 1 


Negative regulation of translation, 
response to stress, negative 
regulation of cell proliferation, 
regulation of elF2 alpha 
phosphorylation by heme 


Enlarged heart size, abnormal 
red blood cell development, 
morphology, physiology with 
macrocytic anemia. 


ns (Ij 


OivrS 


Rnf6 


ring finger protein 
(C3H2C3 type) 6 


Ubiquitin-dependent protein 
catabolic process, positive 
regulation of transcription, 
DNA-dependent 


NA 


ns (2) 


Qivr 16 


Robol 


roundabout homolog 1 
(Drosophila) 


Cell differentiation, axon guidance, 
chemotaxis 


Neonatal death, aphagia, delayed 
lung maturation and bronchial 
hyperplasia. 


insertion 


Qivr 16 


Nrip 1 


Nuclear receptor 
interaacting protein 1 


Regulation of transcription 


Female infertility due to ovulation 
failure. Male and female mice are 
smaller than wild-type littermates. 


ns (1) 


Qivr 16 


Usp2S 


ubiquitin specific 
peptidase 25 


Ubiquitin-dependent protein 
catabolic process 


NA 


ns (1) 


Qivr 17-1 


Cryaa 


crystallin, alpha A 


Negative regulation of apoptosis, 
negative regulation of caspase 
activity, lens fiber cell 
morphogenesis 


Small lenses that develop 
progressive opacity beginning 
in the nucleus. 


no 


Qivr 17-2 


Snfllk/Sikl 


SNFl-like kinase 


Negative regulation of transcription 
from RNA polymerase II promoter, 
regulation of cell differentiation, 
protein kinase cascade 


NA 


no 


Qivr 17-2 


March2 


membrane-associated 
ring finger (C3HC4) 2 


Endocytosis, biological process 


NA 


no 


Qivr 17-2 


Tapbp 


TAP binding protein 


Antigen processing and 
presentation of exogenous peptide 
antigen via MHC class 1, 
TAP-dependent; defense response 


Reduced and thermolabile MHC 
class 1 surface expression due to 
impaired peptide loading with 
stabilizing peptides, impaired 
T cell selection, altered NK 
repertoire, lower CD8+ T cell 
numbers, impaired responses to 
select class l-restricted antigens. 


ns (1) 


Qivr 17-2 


H2-00 


histocompatibility 2, 
0 region alpha locus 


Antigen processing and presentation 
of peptide or polysaccharide antigen 
via MHC class II, regulation of T cell 
differentiation, Graft-versus-host 
disease, viral myocarditis 


Abnormal antigen presentation 
via MHC class II, enhanced 
selection of CD4-h single positive 
thymocytes. Mice homozygous 
for a different knock-out allele 


ns (1) 



show increased serum IgGl levels. 
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Table 2 GO-terms and functions observed in knock-out mice of 31 potential candidates from QTL intervals (Continued) 



Qivrl7-2 H2-DMa histocompatibility 2, Antigen processing and presentation Impaired antigen presenting cell ns (2) 

class II, locus DMa of exogenous peptide antigen via function, poor IgG responses to 









MHC class II, positive regulation of 
T cell differentiation, positive 
regulation of immune response, 
Graft-versus-host disease, viral 
myocarditis 


T-dependent antigens, reduced 
numbers of mature CD4+ T cells, 
increased susceptibility to 
Leishmania major infection. 




Qivrl7-2 


Tap2 


transporter 2, ATP-binding 
cassette, sub-family B 
(MDR^AP) 


Antigen processing and presentation 
of exogenous protein antigen via 
MHC class lb, TAP-dependent; 
positive regulation of T cell mediated 
cytotoxicity, protection from natural 
killer cell mediated cytotoxicity 


No CD8+ T cells, although 
numbers of CD4-I- T cells and 
B cells are normal. 


ns (6) 


Oivrl7-2 


H2-0b 


histocompatibility 2, 
0 region beta locus 


Antigen processing and presentation 
of peptide or polysaccharide antigen 
via MHC class II, Graft-versus-host 
disease, viral myocarditis 


NA 


ns (7) 


Oivrl7-2 


H2-Abl 


histocompatibility 2, 
class II antigen A, beta 1 


Antigen processing and presentation 
of peptide or polysaccharide antigen 
via MHC class II, Graft-versus-host 
disease, viral myocarditis 


Depletion of mature CD4+ T cells, 
deficiency in cell-mediated 
immune responses, increased 
susceptibility to viral infections. 


ns (9) 
stop_L 


Qivrl7-2 


H2-Aa 


histocompatibility 2, 
class II antigen A, alpha 


Antigen processing and presentation 
of exogenous peptide antigen via 
MHC class II, positive regulation of 
T cell differentiation, Graft-versus-host 
disease, viral myocarditis 


Lack of cell surface expression 
of MHC class II molecules on 
macrophages, decreased 
CD4-positive T cell number, 
increased CD8-positive 
T cell number, thymus hyperplasia, 
enlarged lymph nodes, altered 
splenocyte response to 
staphylococcal enterotoxin B. 


ns (10) 


Qivrl7-2 


Lstl 


leukocyte specific 
transcript 1 


Negative regulation of lymphocyte 
proliferation, immune response, 
cell morphogenesis 


NA 


ns (1) 
deletion 


Oivrl7-2 


Gtf2h4 


general transcription 
factor II H, polypeptide 4 


Regulation of transcription, 
DNA-dependent 


NA 


ns (3) 


Oivrl7-2 


H2-T23 


histocompatibility 2, 
T region locus 23 


Antigen processing and presentation 
of peptide antigen via MHC class 1, 
Graft-versus-host disease, viral 
myocarditis 


CD4^ T cells have enhanced 
responses after infection or 
immunization, are resistant to 
suppressor activity mediated by 
a subset of CD8^ T cells, but are 
more susceptible to NK cell lysis. 


ns (3) 


Oivrl7-2 


H2-BI 


histocompatibility 2, 
blastocyst 


Antigen processing and presentation 


NA 


ns (4) 
stop_L 


Oivrl7-2 


Rpp21 


ribonuclease P 21 
subunit (human) 


tRNA processing 


NA 


ns (1) 
stop_G 


Oivrl7-2 


Trim26 


tripartite motif 
protein 26 


Biological process 


NA 


ns (1) 


Oivrl7-2 


Plo2g7 


phospholipase A2, 
group VII (platelet-activating 
factor acetyl hydrolase, 
plasma) 


Inflammatory response, lipid catabolic 
process 


NA 


ns (1) 


Qivrl7-2 


Cyp39ol 


cytochrome P450, 
family 39, subfamily a, 
polypeptide 1 


Lipid metabolic process, oxidation 
reduction 


NA 


ns (2) 


Oivr-19 


Sorbs! 


sorbin and SH3 domain 
containing 1 


Transport, focal adhesion assembly, 
positive regulation of establishment 
of protein localization in plasma 
membrane 


Decreased triglyceride levels, 
altered glucose homeostasis, 
decreased white blood cells and 
resistance to developing glucose 
intolerance induced by a high 
fat diet. 


ns (6) 


Oivr-19 


TctnS 


tectonic family member 3 


Apoptosis 


NA 


ns (1) 



insertion (2) 



Nedeiko et al. BMC Genomics 2012, 13:41 1 
http://www.bionnedcentral.conn/1471 -21 64/1 3/41 1 



Page 1 5 of 1 9 



Table 2 GO-terms and functions observed in knock-out mice of 31 potential candidates from QTL intervals (Continued) 

Oivr-19 Hps! Hermansky-Pudlak Positive regulation of natural killer Hypopigmentation and increased ns (6) 

syndrome 1 homolog cell activation, secretion of lysosomal bleeding time. Impaired natural 

(human) enzymes killer cell function, reduced 

secretion of kidney lysosomal 
enzymes, abnormal retinofugal 
neuronal projections 
characterize some alleles. 

Qivr-19 Dnmbp dynamin binding protein Intracellular signaling cascade, NA ns (3) 

regulation of Rho protein signal 
transduction 

The genotype of C57BL/6J was used as reference for the sequence polymorphisms. The knockout mutant phenotype was identified in the MGI Mouse Genome 
Database. Gene Fign (fidgetin) fulfilled most criteria but was found to be expressed only at late times p.i. (after day 14) and was therefore omitted from the 
candidate gene list. Abbreviations: KO: knockout mutant, NA: not analyzed, SNP: single nucleotide polymorphism, ns: non-synonymous codons, ORF: open reading 
frame, stop_G: stop codon gained, stop_L: stop codon lost. 



humans, increased activities of certain variants of PLA2G7 
were associated with early coronary atherosclerosis and 
with endothelial dysfunction, but the gene may also exert 
an anti-inflammatory function [56-60]. The Pla2g7 gene 
was also identified as a potential candidate gene for suscep- 
tibility against infections with H3N2 influenza virus [33]. 
Pla2g7 expression levels in susceptible A/J mice were 
higher than in resistant C57BL/6J mice after infection with 
H3N2 virus [33]. We also showed previously that Pla2g7 
exhibits a cis-eQTL between C57BL/6J and DBA/2J in 
non-infected lungs where the DBA/2J allele shows high 
levels of expression [26]. Tnfrsf21 which was identified by 
[33] as potential candidate of Qivrl7-2 also exhibits a cis- 
eQTL in non-infected BXD mice [26] but was not found to 
be regulated in C57BL/6 mice after infection (data not pub- 
lished). Tapbp (TAP binding protein) plays a major role in 
the antigen processing and MHC class I presentation by 
stabilizing the TAP peptide transporter, e,g, [61-65]. Also, 
Tap2 (transporter 2, ATP-binding cassette, sub-family B 
(MDR/TAP)) gene is involved in antigen processing and 
presentation [63,66]. Gtf2h4 (general transcription factor II 
H, polypeptide 4) encodes a general transcription factor. 
Recruitment and activation of Gtj2h4 represents a rate- 
limiting step for the emergence of HIV from latency and se- 
quence variants have been associated with multiple scler- 
osis [67-69]. 

Within the Qivrl9 interval, only one gene, Hpsl 
(Hermansky-Pudlak syndrome 1 homolog (human)), 
has been associated with the host responses to infection. 
Mice carrying a natural mutation in the Hpsl gene 
showed an increased inflammatory response in alveolar 
macrophages after intranasal challenge with LPS [70] . 

The GeneNetwork database allows searching for other 
phenotypic traits that exhibit a genome-wide significant 
(LRS > 18) within the Qivr intervals identified by our study. 
Two phenotypic traits, related to neuronal responses (trait 
ID 11285) and body weight changes (trait ID 12838), are 
located to the Qivr 16 locus. Also, the Qivr 17-2 interval 
contained significant QTLs for other traits. Two traits are 
related to host infectious diseases, 'Ectromelia virus survival' 



(ID 12672) and 'Chlamydia psittaci (6BC) infection re- 
sponse' (ID 11025) and four traits are associated with im- 
mune cell responses (ID 10201, 10466, 10238, 10236). In 
addition two traits described seizure responses (ID 10388, 
10507), and one trait has not been disclosed yet (ID 13920). 
Within the early time chromosome 10 interval, three other 
traits exhibit their most significant QTLs: '3a,5a-THDOC 
in blood plasma 3 days after cycle 5 of chronic intermittent 
air vapor' (ID 13027) and two non-disclosed traits. The first 
trait may relate to stress responses in the central nervous 
system (ID 13292 and 13846). 

Conclusions 

The mapping of resistance and susceptibility loci in the 
BXD population revealed several new QTLs and poten- 
tial gene candidates that may be critical for the host 
defense against influenza A virus infection. Body weight 
and survival QTLs showed a time-dependent profile in- 
dicating that the genetic factors in these QTLs are im- 
portant for the host response in a temporal dynamic 
fashion. Five QTL regions were examined in detail, and 
we identified several possible candidate genes that may 
be critical for the host response to influenza A infections 
in humans. 

Methods 

Animals 

The mouse inbred strains C57BL/6J, DBA/2J and 
B6D2F1 were delivered from Janvier, France. Recombin- 
ant inbred mouse strains BXD were purchased from three 
different sources: The Jackson Laboratory, the University 
of Tennessee Health Science Center (Memphis, TN) and 
from Harlan, The Netherlands. For the analysis, mice 
were transferred to the animal facility in Braun- 
schweig and adapted for at least two weeks to the 
new environment before starting experiments. Ani- 
mals were maintained under specific pathogen free 
conditions. AU experiments in mice were approved by 
an external committee and according to the national 
guidelines of the animal welfare law in Germany 
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('Tierschutzgesetz in der Fassung der Bekanntma- 
chung vom 18. Mai 2006 (BGBL I S. 1206, 1313), das 
zuletzt durch Artikel 20 des Gesetzes vom 9. Dezem- 
ber 2010 (BGBL I S. 1934) geandert worden ist/). The 
protocol used in these experiments has been reviewed 
by an ethics committee and approved by the 'Nieder- 
sachsisches Landesamt fiir Verbraucherschutz und 
Lebensmittelsicherheit, Oldenburg, Germany', accord- 
ing to the German animal welfare law (Permit Num- 
bers: 33.42502/04-108/06, 33.9.42502-04-051/09). 

Virus and infection of mice 

The mouse-adapted influenza strain A/Puerto Rico/8/34 
(HlNl; PR8M, [18] and references therein) was used 
for all infection studies. Stocks were prepared by infec- 
tion of 10-day-old embryonated chicken eggs. After 
mice were anesthetized by intra-peritoneal injection of 
Ketamin-Xylazine solution in sterile NaCl (50 mg/ml 
Ketamine, Invesa Arzneimittel GmbH, Freiburg; 2% 
Xylazine, Bayer Health-Care, Leverkusen) with a dose 
adjusted to the individual body weight, mice were 
infected intranasally with 2 x 10^ FFU of PR8M in 20 [d 
of sterile phosphate-buffered saline. Mice were assayed 
daily for body weight (determined as % of initial weight 
at day 0) and mortality during 13 days p.i. We used 
death as the end point for survival. Mice were sacrificed 
if body weight loss exceeded 25%. It should be noted 
that for mice that did not show any signs of body weight 
loss over the entire time period after infection, we do 
not have additional parameters to verify that they have 
indeed been infected. However, these cases were very 
few. In addition, we have ample experience with this 
infection method and the failure rate, for example with 
the DBA/2J strain, is less that 5%. 

Data handling and statistical analysis 

In total, 283 BXD mice and 127 mice from the parental 
strains or Fl generation were used for the infection experi- 
ments. In total 53 BXD strains were infected with an aver- 
age of 5 mice per strain. We performed all primary 
calculations using simple features and ftinctions of Micro- 
soft Excel. Three sets of analysis were performed for the fol- 
lowing variables: (1) body weight in percentage from 
starting weight (day 0) using the strain medians to exclude 
outliers; (2) survival by calculating the survival rate of each 
strain from day 7 to 13, (3) mean time to death in days. For 
statistical analyses, tests and visualization we used R, a free 
software environment for statistical computing and graph- 
ics (http://www.R-project.org). In order to test for batch 
effects or other co-factors, we visualized the data using 
multidimensional scaling based on the Sammon mapping 
method [71]. No clusters with respect to any of the co- 
factors age, weight at day 0, experimenter, time of infection, 
or source of mice could be found in the visualizations. 



QTL mapping 

QTL mapping was performed using the web-based com- 
plex trait analysis available on the GeneNetwork website 
(www.genenetwork.org) and the mapping module to 
analyze phenotypes in context of mouse genotypic differ- 
ences. Interval mapping evaluates potential QTLs at 
regular intervals and estimates the significance at each 
location with a graphical representation of the likelihood 
ratio statistics (LRS) using 2000 permutation tests 
[19,22]. LRS values may be converted to LOD scores by 
dividing by 4.61. For the two locus model the following 
equation was used: Var =Q1 -\-Q2-\- QlxQ2 + e, where 
Var = the between-strain mean variance in the trait, Ql 
and Q2 are makers tightly linked to the loci, Q1XQ2 
is the additive-by-additive' epistatic interaction term, 
and e is the residual error. The original data sets can be 
obtained at www.genenetwork.org with the following 
identification numbers: body weight: 13005 to 13017; 
survival: 13000 to 13004 and mean time to death: 12996. 
We performed full genome scans for epistatic interac- 
tions using the pair-scan module that is implemented in 
GeneNetwork. This module exploits the direct global 
optimization algorithm developed by [72]. The code 
compares the fit (as measured by LRS scores) for a 
purely additive model, a purely epistatic model, and the 
full model. The code also implements a permutation test 
(n = 500) and this enabled us to estimate the empirical p 
value of the alternative models. 

Candidate gene discovery 

The QTL region analysis was initially performed 
using the QTLminer which has been implemented in 
GeneNetwork [25]. By using the automatic function 
of GeneNetwork we identified significant cis-QTLs 
with LRS higher than 18 at a genome-wide p-value 
of < 0.05. Additionally the genes mapped within the 
analyzed QTLs were surveyed by the National Center 
for Biotechnology Information (NCBI) Entrez Gene 
website (http://www.ncbi.nlm.nih.gov/sites/entrez?db= 
gene) and the Jackson Laboratory's MGI Mouse 
Genome Database project (http://www.informatics. 
jax.org/) to identify potential candidate genes. The 
GeneRIF database (http://www.ncbi.nlm.nih.gov/pro- 
jects/GeneRIF/GeneRIFhelp.html) was used as a pri- 
mary source to search for known gene functions and 
corresponding citations. 

All sequence variants between B6 and D2 parental 
genomes (SNPs, indels) were extracted by using a 
comparative analysis that relies on approximately lOOx 
whole genome shot gun of DBA/2J [73]. All of these 
sequence data are available at http://ucscbrowser. 
genenetwork.org/, the GeneNetwork Variant Browser 
(http://www.genenetwork.org/webqtl/main.py?FormID= 
snpBrowser), and the NCBI Short Read Archive (18 files 
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total, e.g., SRX037575, SRX013980, SRX013299, SRX012582, 
SRX012581, SRX012580); http://www.biomedcentraLcom/ 
1471-2105/11/S4/O7. 



Additional files 



Additional file 1: Table SI. Body weight loss from day 1 until day 7 for 
all strains showing the mean values, SEM per day and number of mice 
per strain analyzed. 

Additional file 2: Figure SI. Differential susceptibility to influenza A 
infection among different BXD strains, parental strains and Fl generation. 
Rank-ordered strain distribution pattern illustrating the percentage of 
surviving mice per strain after influenza A infection for 53 BXD, parental 
DBA/2J and C57BL/6J strains and Fl (B6D2F1) mice for day 7 (A, trait ID: 
13001) and day 8 (B, trait ID: 13002). Rank-ordered strain distribution 
pattern showing mean time to death for 53 BXD, parental DBA/2J and 
C57BL/6J strains and Fl (B6D2F1) mice (C, trait ID: 12996). 

Additional file 3: Figure S2. Principal component analysis of all body 
weight traits. Percent contribution of principal components to the total 
variance for body weight traits from day 1 until day 7 p.i. PCI contributes 
about 80% and PC2 about 10% to the total variance. 

Additional file 4: Figure S3. Body weight changes in mock-infected 
C57BL/6J and DBA/2J mice. Female DBA/2J (n=4) and C57BL/6J (n=3) 
mice were intranasally inoculated with 25 [i\ PBS under anesthesia. Body 
weight changes for each group of treated mice at various days p.i. is 
shown with reference to the starting weight (% body weight). Data 
represent mean values +/- SEM. 

Additional file 5: Figure S4. Genome-wide linkage analysis for survival 
after HlNl infection at days 9 and 10 p.i. Analysis of survival from day 9 
(A) and day 10 (B) after PR8 infection revealed suggestive QTLs on 
chromosomes 5, 16, 17 and 19. Interval mapping is shown across the 
entire genome as blue line representing the LRS scores; green line: DBA/ 
2J alleles increase trait values; red line: C57BL/6J alleles increase trait 
values. Upper x-axis shows mouse chromosomes, lower x-axis shows 
physical map in mega bases for each chromosome, y-axis represents LRS 
score. The horizontal lines mark the genome-wide significant thresholds 
at p0.05 (red line) and suggestive thresholds at pO.37 (gray line). 

Additional file 6: Figure S5. Pair Scan analysis of trait 13002. The table 
provides the breakdown of trait 13002 as an example for a pair scan 
analysis study that compares the fit of the alternative models. An analysis 
of this type can be regenerated rapidly on any of the traits presented in 
our manuscript. 

Additional file 7: Figure S6. Scheme to identify potential candidate 
genes in the analyzed QTL intervals. The detailed steps are described in 
the main manuscript. 

Additional file 8: Table S2. List of genes in the QTL intervals that were 
studied in detail. The QTLminer tool in GeneNetwork was to identif/ all 
genes in the QTL intervals and to obtain associated GO terms and KEGG 
pathways. All annotated genes within the QTL intervals that were 
expressed between day 1 and day 60 after infection of C57BL/6J mice 
with PR8M virus were further analyzed for the following attributes: genes 
that were up- or down-regulated in infected lungs after PR8M infection, 
genes carrying an insertion or deletion or a non-synonymous nucleotide 
change in the open reading frame, genes exhibiting a cis-expression QTL 
in the lung, and genes that were differentially expressed in C57BL/6J and 
DBA/2J mice in infected lungs between day 1 and day 8 p.i. 
Abbreviations: FS: frame shift in coding region, SNP: single nucleotide 
polymorphism, cis-eQTL: cis-regulated expression QTL, B6: C57BL/6J, D2: 
DBA/2J, ns: non-synonymous codons, stop_G: stop codon gained, stop_L: 
stop codon lost, DEL: deletion, INS: insertion; ns (. . .): number of non- 
synonymous codons. Legends: Expression d0-d60, log2 >8 at any day 
during infection; No: log2 8; ?: no expression data available. Regulation 
during infection (d0-d60): + - up or down regulation, genes that 
exhibited at least a 2-fold difference in changes of the expression levels 
in lungs of infected mice compared to non-infected controls; (+)(-): up-or 
down-regulation, genes that exhibited at least a 



1.5-fold difference in changes of the expression levels in lungs of infected 
mice compared to the non-infected controls; no: no regulation. SNP: ns 
(..) number of non-synonymous codons; STOP_L: Stop codon lost, 
STOP_G: Stop codon gained; no: no SNPs between B6 and D2. Indel: DEL: 
Deletion; INS: Insertion; no: no INDEL. cisQTL: yes (. . .): LRS value has to 
be greater than >18; no: no QTL or LRS 18. diff expressed (B6-D2) during 
infection: Genes were defined as differentially expressed (log2 >8) at any 
day in infected lungs of C57BL/6J and DBA/2J; + -: Genes were defined 
as differentially expressed exhibiting at least a 2.0 difference in changes 
of the expression levels in infected lungs of C57BL/6J and DBA/2J; (+)(-) 
differentially expressed: Genes were defined as differentially expressed 
exhibiting at least a 1.5 difference in changes of the expression levels in 
infected lungs of C57BL/6J and DBA/2J; no: no regulation or not 
expressed; Score: the different attributes were counted. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

TN, MH and KS conceived and designed the experiments. TN and MH 
infected the mice and collected the data. HK, RA, RW and KS performed and 
interpreted the QTL analysis. FK performed the analysis of environmental and 
experimental factors. SS and LL provided BXD mice for experimental studies. 
HK, RW and KS wrote the manuscript. All authors read and approved the 
final manuscript. 

Authors' information 

Tatiana Nedeiko and Heike Kollmus contributed equally as first authors. 
Robert W. Williams and Klaus Schughart contributed equally as senior 
authors. 

Acknowledgments 

This work was supported by intra-mural grants from the Helmholtz- 
Association (Program Infection and Immunity) and a research grant 
FluResearchNet (No. 01KI07137) from the German Ministry of Education and 
Research to KS and the virtual institute 'GeNeSys' funded by the Helmholtz 
Association. The funders have no role in study design, data collection and 
analysis, decision to publish, or preparation of the manuscript Mice for these 
experiments were maintained by the animal caretakers of the Central Animal 
Facilities at the HZI. We wish to thank Christin Fricke for excellent technical 
assistance. 

Author details 

^Department of Infection Genetics, Helmholtz Centre for Infection Research 
and University of Veterinary Medicine Hannover, 38124, Braunschweig, 
Germany. ^Department of Bioinformatics and Statistics, Helmholtz Centre for 
Infection Research, Braunschweig, Germany. ^Department of Computer 
Science, Ostfalia University of Applied Sciences, Wolfenbuttel, Germany. 
^Department of Anatomy and Neurobiology, University of Tennessee Health 
Science Center, Memphis, Tennessee, United States of America. ^Department 
of Molecular and Cellular Neurobiology, Neuroscience Campus Amsterdam, 
Amsterdam, VU, the Netherlands. ^Jiangsu Key Laboratory of 
Neuroregeneration, Nantong University, Nantong, China. ^Nycomed GmbH, 
Institute for Pharmacology and Preclinical Drug Safety, 
Barsbuettel-Willinghusen, Germany. 

Received: 10 April 2012 Accepted: 10 August 2012 
Published: 20 August 2012 

References 

1. Johnson NP, Mueller J: Updating the accounts: global mortality of the 
1918-1920 "Spanish" influenza pandemic. Bull Hist Med 2002, 
76(1):105-115. 

2. Fauci AS: Seasonal and pandemic influenza preparedness: science and 
countermeasures. J Infect Dis 2006, 194(Suppl 2):S73-76. 

3. Scriven J, McEwen R, Mistry S, Green C, Osman H, Bailey M, Ellis C: 
Swine flu: a Birmingham experience, Clin Med 2009, 9(6):534-538. 

4. Yates L, Pierce M, Stephens S, Mill A, Spark P, Kurinczuk J, Valappil M, 
Brocklehurst P, Thomas S, Knight M: Influenza A/HlNlv in pregnancy: an 
investigation of the characteristics and management of affected women 



Nedeiko et al. BMC Genomics 2012, 13:41 1 
http://www.bionnedcentral.conn/1471 -21 64/1 3/41 1 



Page 18 of 19 



and the relationship to pregnancy outcomes for mother and infant. 

Health Technol Assess 2010, 14(34):! 09-1 82. 

5. Albright FS, Orlando P, Pavia AT, Jackson GG, Cannon-Albright LA: Evidence 
for a heritable predisposition to death due to influenza. J Infect Dis 2008, 
1 97(1 ):1 8-24. 

6. Gottfredsson M, Halldorsson BV, Jonsson S, Kristjansson M, Kristjansson K, 
Kristinsson KG, Love A, Blondal T, Viboud C, Thorvaldsson S, et al: Lessons 
from the past: familial aggregation analysis of fatal pandemic 
influenza (Spanish flu) in Iceland in 1918. Proc Natl Acad Sci USA 2008, 
105(4):1 303-1 308. 

7. Horby P, Sudoyo H, Viprakasit V, Fox A, Thai PQ, Yu H, Davila S, Hibberd M, 
Dunstan SJ, Monteerarat Y, et al: What is the evidence of a role for host 
genetics in susceptibility to influenza A/H5N1? Epidemiol Infect 2010, 
138(1 l):l-9. 

8. Everitt AR, Clare S, Pertel T, John SP, Wash RS, Smith SE, Chin CR, Feeley EM, 
Sims JS, Adams DJ, et al: IFITM3 restricts the morbidity and mortality 
associated with influenza. Nature 2012, 484(7395):5 19-523. 

9. Srivastava B, Blazejewska P, Hessmann M, Bruder D, Geffers R, Mauel S, 
Gruber AD, Schughart K: Host genetic background strongly influences the 
response to influenza a virus infections. PLoS ONE 2009, 4(3):e4857. 

10. Trammell RA, Toth LA: Genetic susceptibility and resistance to influenza 
infection and disease in humans and mice. Expert Rev Mol Diagn 2008, 
8(4):5 15-529. 

11. Ding M, Lu L, Toth LA: Gene expression in lung and basal forebrain 
during influenza infection in mice. Genes Brain Behav 2008, 7(2):173-183. 

12. Boon AC, DeBeauchamp J, Hollmann A, Luke J, Kotb M, Rowe S, 
Finkelstein D, Neale G, Lu L, Williams RW, et al: Host genetic variation 
affects resistance to infection with a highly pathogenic H5N1 influenza 
A virus in mice. J Virol 2009, 83(20):1 041 7-1 0426. 

13. Boon AC, Debeauchamp J, Krauss S, Rubrum A, Webb AD, Webster RG, 
McElhaney J, Webby RJ: Cross-reactive neutralizing antibodies directed 
against pandemic H1N1 2009 virus are protective in a highly sensitive 
DBA/2 influenza mouse model. J Virol 2010, 84(15):7662-7667. 

14. Otte A, Sauter M, Alleva L, Baumgarte S, Klingel K, Gabriel G: Differential 
host determinants contribute to the pathogenesis of 2009 pandemic 
H1N1 and human H5N1 influenza A viruses in experimental mouse 
models. Am J Pathol 201 1, 179(l):230-239. 

15. Boon AC, Finkelstein D, Zheng M, Liao G, Allard J, Klumpp K, Webster R, 
Peltz G, Webby RJ: H5N1 influenza virus pathogenesis in genetically 
diverse mice is mediated at the level of viral load. MBio 201 1, 2(5):pii: 
e00171-ll. 

16. Trammell RA: Liberati TA. Toth LA: Host genetic background and the innate 
inflammatory response of lung to influenza virus. Microbes Infect; 201 1. 

17. Pica N, Iyer A, Ramos I, Bouvier NM, Fernandez-Sesma A, Garcia-Sastre A, 
Lowen AC, Palese P, Steel J: The DBA.2 mouse is susceptible to disease 
following infection with a broad, but limited, range of influenza A and B 
viruses. J Virol 201 1, 85(23):5-12829. 

18. Blazejewska P, Koscinski L, Viegas N, Anhlan D, Ludwig S, Schughart K: 
Pathogenicity of different PR8 influenza A virus variants in mice 

is determined by both viral and host factors. Virology 201 1, 
412(l):36-45. 

1 9. Peirce JL, Lu L, Gu J, Silver LM, Williams RW: A new set of BXD 
recombinant inbred lines from advanced intercross populations in mice. 
BMC Genet 2004, 5:7. 

20. Taylor BA, Wnek C, Kotlus BS, Roeme RN, MacTaggart T, Phillips SJ: 
Genotyping new BXD recombinant inbred mouse strains and 
comparison of BXD and consensus maps. Mamm Genome 1999, 
10(4):335-348. 

21. Williams RW, Gu J, Qi S, Lu L: The genetic structure of recombinant inbred 
mice: high-resolution consensus maps for complex trait analysis. Genome 
Biol 2001, 2(11):RESEARCH0046. 

22. Chesler EJ, Lu L, Wang J, Williams RW, Manly KF: WebQTL: rapid 
exploratory analysis of gene expression and genetic networks for brain 
and behavior. Nat Neurosci 2004, 7(5):485-486. 

23. Alberts R, Srivastava B, Wu H, Viegas N, Geffers R, Klawonn F, Novoselova N, 
Do Valle TZ, Panthier JJ, Schughart K: Gene expression changes in the 
host response between resistant and susceptible inbred mouse strains 
after influenza A infection. Microbes Infect 2010, 12(4):309-318. 

24. Ziebarth JD, Cook MN, Wang X, Williams RW, Lu L, Cui Y: Treatment- and 
population-dependent activity patterns of behavioral and expression 



QTLs. PLoS ONE 2012, 7(2):e31805. 

25. Alberts R, Schughart K: QTLminer: identifying genes regulating 
quantitative traits. BMC Bioinformatics 2010, 1 1:516. 

26. Alberts R, Lu L, Williams RW, Schughart K: Genome-wide analysis of the 
mouse lung transcriptome reveals novel molecular gene interaction 
networks and cell-specific expression signatures. Respir Res 201 1, 12:61. 

27. Zhirnov OP, Matrosovich TY, Matrosovich MN, Klenk HD: Aprotinin, a 
protease inhibitor, suppresses proteolytic activation of pandemic HINIv 
influenza virus. Antivir Chem Chemother 201 1 , 21 (4):1 69-1 74. 

28. Okumura Y, Takahashi E, Yano M, Ohuchi M, Daidoji T, Nakaya T, Bottcher E, 
Garten W, Klenk HD, Kido H: Novel type II transmembrane serine 
proteases, MSPL and TI\/1PRSS13, proteolytically activate membrane 
fusion activity of hemagglutinin of highly pathogenic avian 
influenza viruses and induce their multicycle replication. J Viro 2010, 
84(10):5089-5096. 

29. Bottcher-Friebertshauser E, Freuer C, Sielaff F, Schmidt S, Eickmann M, 
Uhlendorff J, Steinmetzer T, Klenk HD, Garten W: Cleavage of influenza 
virus hemagglutinin by airway proteases TI\/1PRSS2 and HAT differs in 
subcellular localization and susceptibility to protease inhibitors. J Virol 

2010, 84(11):5605-5614. 

30. Bertram S, Glowacka I, Blazejewska P, Soilleux E, Allen P, Danisch S, Steffen I, 
Choi SY, Park Y, Schneider H, et al: TI\/1PRSS2 and TI\/1PRSS4 facilitate 
trypsin-independent spread of influenza virus in Caco-2 cells. J Virol 2010, 
84(1 9):1 001 6-1 0025. 

31. Chaipan C, Kobasa D, Bertram S, Glowacka I, Steffen I, Tsegaye TS, Takeda M, 
Bugge TH, Kim S, Park Y, et al: Proteolytic activation of the 1918 influenza 
virus hemagglutinin. J Virol 2009, 83(7):3200-321 1. 

32. Steinhauer DA: Role of hemagglutinin cleavage for the pathogenicity of 
influenza virus. Virology 1999, 258(1 ):1 -20. 

33. Boivin GA, Pothlichet J, Skamene E, Brown EG, Loredo-Osti JC, Sladek R, 
Vidal SM: Mapping of clinical and expression quantitative trait loci in a 
sex-dependent effect of host susceptibility to mouse-adapted influenza 
H3N2/HK/1/68. J Immunol 2012, . 

34. Bottomly D, Ferris MT, Aicher LD, Rosenzweig E, Whitmore A, Aylor DL, 
Haagmans BL, Gralinski LE, Bradel-Tretheway BG, Bryan JT, et al: Expression 
quantitative trait Loci for extreme host response to influenza a in 
pre-collaborative cross mice. G3 (Bethesda) 2012, 2(2):21 3-221. 

35. Krishnamoorthy J, Mounir Z, Raven JF, Koromilas AE: The elF2alpha kinases 
inhibit vesicular stomatitis virus replication independently of elF2alpha 
phosphorylation. Cell Cycle 2008, 7(15):2346-2351. 

36. Thakur SA, Zaiinger Z, Johnson TR, Imani F: PKR is a novel mediator 
of CD40 signaling and plays a critical role in modulating 
immunoglobulin expression during RSV infection. Clin Vaccine Immunol 

2011, 18(12):2060-2066. 

37. Bergmann M, Garcia-Sastre A, Carnero E, Pehamberger H, Wolff K, Palese P, 
Muster T: Influenza virus NS1 protein counteracts PKR-mediated 
inhibition of replication. J Virol 2000, 74(13):6203-6206. 

38. Ludlow A, Yee KO, Lipman R, Bronson R, Weinreb P, Huang X, Sheppard D, 
Lawler J: Characterization of integrin beta6 and thrombospondin-1 
double-null mice. J Cell Mol Med 2005, 9(2):421-437. 

39. Hogmalm A, Sheppard D, Lappalainen U, Bry K: beta6 Integrin subunit 
deficiency alleviates lung injury in a mouse model of bronchopulmonary 
dysplasia. Am J Respir Cell Mol Biol 2010, 43(l):88-98. 

40. Yoneyama M, Kikuchi M, Matsumoto K, Imaizumi T, Miyagishi M, Taira K, 
Foy E, Loo YM, Gale M Jr, Akira S, et al: Shared and unique functions of 
the DExD/H-box helicases RIG-I, MDA5, and LGP2 in antiviral innate 
immunity. J Immunol 2005, 175(5):285 1-2858. 

41. Choi MK, Wang Z, Ban T, Yanai H, Lu Y, Koshiba R, Nakaima Y, Hangai S, 
Savitsky D, Nakasato M, et al: A selective contribution of the RIG-l-like 
receptor pathway to type I interferon responses activated by cytosolic 
DNA. Proc Natl Acad Sci USA 2009, 1 06(42):1 7870-1 7875. 

42. Gitlin L, Benoit L, Song C, Cella M, Gilfillan S, Holtzman MJ, Colonna M: 
Melanoma differentiation-associated gene 5 (I\/1DA5) is involved in the 
innate immune response to Paramyxoviridae infection in vivo. PLoS 
Pathog 2010, 6(1 ):el 000734. 

43. Gitlin L, Barchet W, Gilfillan S, Cella M, Beutler B, Flavell RA, Diamond MS, 
Colonna M: Essential role of mda-5 in type I IFN responses to 
polyriboinosinic:polyribocytidylic acid and encephalomyocarditis 
picornavirus. Proc Natl Acad Sci USA 2006, 1 03(22):8459-8464. 

44. Kato H, Takeuchi 0, Sato S, Yoneyama M, Yamamoto M, Matsui K, 



Nedeiko et al. BMC Genomics 201 2, 1 3:41 1 Page 1 9 of 1 9 

http://www.bionnedcentral.conn/1471 -21 64/1 3/41 1 



Uematsu S, Jung A, Kawai T, Ishii KJ, et al: Differential roles of MDA5 
and RIG-I helicases in the recognition of RNA viruses. Nature 2006, 
441 (7089):! 01 -105. 

45. Andrews W, Barber M, Hernadez-Miranda LR, Xian J, Rakic S, Sundaresan V, 
Rabbitts TH, Pannell R, Rabbitts P, Thompson H, et al: The role of Slit-Robo 
signaling in the generation, migration and morphological differentiation 
of cortical interneurons. Dev Biol 2008, 313(2):648-658. 

46. Long H, Sabatier C, Ma L, Plump A, Yuan W, Ornitz DM, Tamada A, 
Murakami F, Goodman CS, Tessier-Lavigne M: Conserved roles for Slit and 
Robo proteins in midline commissural axon guidance. Neuron 2004, 
42(2):21 3-223. 

47. Xian J, Clark KJ, Fordham R, Pannell R, Rabbitts TH, Rabbitts PH: Inadequate 
lung development and bronchial hyperplasia in mice with a targeted 
deletion in the Duttl/Robol gene. Proc Natl Acad Sci USA 2001, 
98(26):1 5062-1 5066. 

48. Clark K, Hammond E, Rabbitts P: Temporal and spatial expression of two 
isoforms of the Duttl/Robol gene in mouse development. FEB5 Lett 
2002, 523(1 -3):1 2-1 6. 

49. Zschiedrich I, Hardeland U, Krones-Herzig A, Berriel Diaz M, Vegiopoulos A, 
Muggenburg J, Sombroek D, Hofmann TG, Zawatzky R, Yu X, et al: 
Coactivator function of RIP140 for NFkappaB/RelA-dependent cytokine 
gene expression. Blood 2008, 1 1 2(2):264-276. 

50. Ho PC, Chang KC, Chuang YS, Wei LN: Cholesterol regulation of receptor- 
interacting protein 140 via microRNA-33 in inflammatory cytokine 
production. FA5EB J 20]], 25(5):1 758-1 766. 

51. Janeway CA, Travers P, Walport MJ, Shiomchik MJ: Immunobiology: The 
Immune System in Health and Disease. New York: Garland Science; 2005. 

52. Ovsyannikova IG, Vierkant RA, Pankratz VS, Jacobson RM, Poland GA: 
Extended LTA, TNF, LST1 and HLA gene haplotypes and their association 
with rubella vaccine-induced immunity. PLoS ONE 2010, 5(7):e1 1806. 

53. Schiller C, Nitschke MJ, SeidI A, Kremmer E, Weiss EH: Rat monoclonal 
antibodies specific for LST1 proteins. Hybridoma (Larchmt) 2009, 
28(4):281-286. 

54. Mulcahy H, O'Rourke KP, Adams C, Molloy MG, O'Gara F: LST1 and NCR3 
expression in autoimmune inflammation and in response to IFN-gamma, 
LPS and microbial infection. Immunogenetics 2006, 57(12):893-903. 

55. Kilding R, lies MM, Timms JM, Worthington J, Wilson AG: Additional genetic 
susceptibility for rheumatoid arthritis telomeric of the DRB1 locus. 
Arthritis Rheum 2004, 50(3):763-769. 

56. Lavi S, McConnell JP, Rihal CS, Prasad A, Mathew V, Lerman LO, Lerman A: 
Local production of lipoprotein-associated phospholipase A2 and 
lysophosphatidylcholine in the coronary circulation: association with 
early coronary atherosclerosis and endothelial dysfunction in humans. 
Circulation 2007, 1 15(21 ):271 5-2721. 

57. Kim JY, Hyun YJ, Jang Y, Lee BK, Chae JS, Km SE, Yeo HY, Jeong TS, 
Jeon DW, Lee JH: Lipoprotein-associated phospholipase A2 activity is 
associated with coronary artery disease and markers of oxidative stress: 
a case-control study. Am J Clin Nutr 2008, 88(3):630-637. 

58. Kostopanagiotou G, Avgerinos E, Costopanagiotou C, Arkadopoulos N, 
Andreadou I, Diamantopoulou K, Lekka M, Smyrniotis V, Nakos G: Acute 
lung injury in a rat model of intestinal ischemia-reperfusion: the 
potential time depended role of phospholipases A(2). J Surg Res 2008, 
1 47(1 ):1 08-1 16. 

59. Sutton BS, Crosslin DR, Shah SH, Nelson SC, Bassil A, Hale AB, Haynes C, 
Goldschmidt-Clermont PJ, Vance JM, Seo D, et al: Comprehensive genetic 
analysis of the platelet activating factor acetylhydrolase (PLA2G7) gene 
and cardiovascular disease in case-control and family datasets. Hum Mol 
Genet 2008, 17(9):1318-1328. 

60. Wang Q, Hao Y, Mo X, Wang L, Lu X, Huang J, Cao J, Li H, Gu D: PLA2G7 
gene polymorphisms and coronary heart disease risk: a meta-analysis. 
Thromb Res 2010, 1 26(6):498-503. 

61. Sadegh-Nasseri S, Chen M, Narayan K, Bouvier M: The convergent roles of 
tapasin and HLA-DM in antigen presentation. Trends Immunol 2008, 
29(3):141-147. 

62. Turnquist HR, Petersen JL, Vargas SE, Mcllhaney MM, Bedows E, Mayer WE, 
Grandea AG 3rd, Van Kaer L, Solheim JC: The Ig-like domain of tapasin 
influences intermolecular interactions. J Immunol 2004, 172(5):2976-2984. 

63. Garbi N, Tiwari N, Momburg F, Hammerling GJ: A major role for tapasin as 
a stabilizer of the TAP peptide transporter and consequences for MHC 
class I expression. Eur J Immunol 2003, 33(1):264-273. 



64. Williams AP, Peh CA, Purcell AW, McCluskey J, Elliott T: Optimization of the 
MHC class I peptide cargo is dependent on tapasin. Immunity 2002, 
16(4):509-520. 

65. Raghuraman G, Lapinski PE, Raghavan M: Tapasin interacts with the 
membrane-spanning domains of both TAP subunits and enhances 
the structural stability of TAP1 x TAP2 Complexes. J Biol Chem 2002, 
277(44):41 786-41 794. 

66. Rufer E, Leonhardt RM, Knittler MR: Molecular architecture of the 
TAP-associated MHC class I peptide-loading complex. J Immunol 2007, 
179(9):571 7-5727. 

67. Fukuda A, Tokonabe S, Hamada M, Matsumoto M, Tsukui T, Nogi Y, 
Hisatake K Alleviation of PC4-mediated transcriptional repression by the 
ERCC3 helicase activity of general transcription factor TFIIH. J Biol Chem 
2003, 278(17):14827-14831. 

68. Km YK, Bourgeois CF, Pearson R, Tyagi M, West MJ, Wong J, Wu SY, 
Chiang CM, Karn J: Recruitment of TFIIH to the HIV LTR is a rate-limiting 
step in the emergence of HIV from latency. EMBO J 2006, 
25(15):3596-3604. 

69. Briggs FB, Goldstein BA, McCauley JL, Zuvich RL, De Jager PL, Rioux JD, 
Ivinson AJ, Compston A, Hafler DA, Hauser SL, et al: Variation within DNA 
repair pathway genes and risk of multiple sclerosis. Am J Epidemiol 2010, 
172(2):21 7-224. 

70. Young LR, Borchers MT, Allen HL, Gibbons RS, McCormack FX: Lung- 
restricted macrophage activation in the pearl mouse model of 
Hermansky-Pudlak syndrome. J Immunol 2006, 176(7):436 1-4368. 

71. Sammon JW: A nonlinear mapping for data structure analysis. IEEE 
Transactions on Computers 1969, 18:401-409. 

72. Ljungberg K, Holmgren S, Carlborg 0: Simultaneous search for multiple 
QTL using the global optimization algorithm DIRECT. Bioinformatics 2004, 
20(1 2):1 887-1 895. 

73. Wang X, Agarwala R, Capra J, Chen Z, Church D, Ciobanu D, Li Z, Lu L, 
Mozhui K, Mulligan M, et al: High-throughput sequencing of the DBA/2J 
mouse genome. BMC Bioinformatics 2010, 1 1(Suppl 4):07. doi:10.1 1 86/1 471- 
2105-1 111-S1 184-01 187. 



doi:1 0.1 1 86/1 471 -21 64-1 3-41 1 

Cite this article as: Nedeiko et al.: Distinct gene loci control the host 
response to influenza H1N1 virus infection in a time-dependent manner. 

BMC Genomics 2012 13:411. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



